Edge states and topological orders in the spin liquid phases of star lattice 
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A group of novel materials can be mapped to the star lattice, which exhibits some novel physical 
properties. We give the bulk-edge correspondence theory of the star lattice and study the edge 
states and their topological orders in different spin liquid phases. The bulk and edge-state energy 
structures and Chern number depend on the spin liquid phases and hopping parameters because the 
local spontaneous magnetic flux in the spin liquid phase breaks the time reversal and space inversion 
symmetries. We give the characteristics of bulk and edge energy structures and their corresponding 
Chern numbers in the uniform, nematic and chiral spin liquids. In particular, we obtain analytically 
the phase diagram of the topological orders for the chiral spin liquid states 5'L[(/), 0, —20], where 
(j) is the magnetic flux in two triangles and a dodecagon in the unit cell. Moreover, we find the 
topological invariance for the spin liquid phases, 5'L[(/)i, 02, — (0i + 02)] and *S'L[02, 0i, — (0i + 02)]. 
The results reveal the relationship between the energy-band and edge-state structures and their 
topological orders of the star lattice. 



I. INTRODUCTION 

The discovery of the Integer Quantum Hall effect 
(IQHE) stimulates novel fundamental concepts in con- 
densed matter physics, such as the gauge invariance, 
edge state, ^ and Chern numbei^. In particular, Y. Hat- 
sugai reveals the relationship between Chern number and 
edge states in the IQHE,^ ^ which provides another way 
from the edge state to understand the topological or- 
der in finite systems.^ Moreover, the successful synthe- 
sis of nano and novel materials, such as the triangular 
organic material hi - BEDT{CN)sJ the ka gome lat- 
tice herbertsmithite^ and the three-dimensional hyper- 
kagome lattice magnet {Na^IrsOs)^^ provide many op- 
portunities to examine theoretically and experimentally 
some novel physical properties, including the topolog- 
ical properties fractionalized excitationP^'^, singlet 
valence-bond solid states ,^1^^^, and edge states. The- 
oretically, these materials can be mapped to novel ge- 
ometric lattices, such as the star lattice which is also 
called the triangle-honeycomb latticd^^^. Fisher lattice 
or decorated honeycomb lattice-fl^ These geometric lat- 
tice models provide a new view to understand the ge- 
ometric effect and spin frustration.'^E^E^ In particular, 
the spin models on the star and honeycomb lattices have 
shown many novel phases including the Abelian and 
non-Abelian anyons, chiral spin-liquid phases, topolog- 
ical orders^fl^^ magnetic orders^^^^, and topological 
insulator.™^ 

The IQHE reveals some novel transport properties of 
electrons in two-dimensional (2D) systems, in which the 
edge states play a key role and the Hall conductance 
can be expressed in terms of Chern number.^ ^ Inter- 
estingly, the bulk-edge correspondence discovered by Y. 
HatsugaPI^ becomes an efficient method to explore the 
edge states in various 2D systems, such as the honey- 
comb lattice (Graphene)(^ and spin-chiral ferromagnetic 
kagome lattice.'^^ However, the star lattice consists of a 
special lattice geometry, which can be mapped to a class 




FIG. 1: (Color online) A star lattice with basis ei and 62. 
There are six sites in the unit cell, tc represents the hopping 
amplitude inside the triangle (the bond in black), and td cor- 
responds to the hopping amplitude between different triangles 
(the bond in blue). 



of materials and cold atoms in optical lattices. The mean 
field study of Heisenberg model demonstrates the exis- 
tence of several spin liquid phases, which depend on the 
flux configurations of two triangles and one dodecagon in 
the unit cell.^^ A natural question arises: what is the re- 
lationship between the spin liquid phases and edge states 
on the star lattice with boundaries? 

In this paper, we focus on the edge states and their 
topological orders on the star lattice with boundaries. 
We begin with a 2D tight-binding model with the Hund's 
rule coupling. It can be mapped to an effective tight- 
binding spinless model. The mean field approach pre- 
dicts that there exists several spin liquid phases in the 
ground states with the local time reversal symmetry 
breaking. We use the bulk-edge correspondence method 
to analyze the edge states and their topological orders on 
the star lattice with boundaries. 

This paper is organized as follows. In Sec. II, we intro- 
duce the tight-binding model with boundaries and map 
it to a spinless tight-binding model. In Sec. HI, we give 
the bulk-edge correspondence for the star lattice. We 
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FIG. 2: (Color online) The elementary plaquette of star lat- 
tice contains two inequivalent triangles <, >, and one do- 
decagon. The magnetic flux configurations are labeled by 
6'L[(/)<],0>,(/)i2] following Ref.^^. 

present the edge states and their corresponding Chern 
numbers in various phases in Sec. IV. Finally, we give 
the discussion and conclusions. 
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FIG. 3: (Color online) The intersection number between the 
canonical loop on the complex-energy surface (Riemann sur- 
face) and the trace of the edge state energy Ref.^. 



II. MODEL AND SPIN LIQUID PHASES 



In order to understand the relationship between the 
lattice geometry, edge states, and their topological prop- 
erties, we consider the star lattice with boundaries, in 
which the conducting electrons move in a local spin 
background and couple with them by the Hund's rule 
to form a double-exchange system. The corresponding 
tight-binding Hamiltonian can be written as. 



H 



H.c.) - J ^ 4a^«/3 • SiCi^ (1) 



where tij is the hopping amplitude between two nearest 
neighboring sites (i, j); ct^(Qcr) is the creation (annihila- 
tion) operator on site i with spin a; Si is the local spin on 
site i, which couples with the conducting electron spins 
with the effective coupling constant J. We consider that 
the local spins are approximately classical and the cou- 
pling J is strong enough to have the hoping electrons 
to align them to the local spin Si on each site with the 
spinon function, |x) = (e^^ cos(l9i/2), e^(^^+^^) sin((9i/2)), 
where {Oi^(j)i) are the spinon parameters. In this spinon 
representation, the Hamiltonian Eq.Q can be mapped 
to an effective tight-binding Hamiltonian, 



H.c.) 



(2) 



dodecagon plaquette. The mean field study of Heisen- 
berg model on the star lattice has revealed the existence 
of several spin liquid phases.^^ The spin liquid phases 
depend on the flux figuration of the unit cell, which 
is labeled by the notation ^12].'^ In terms 

of the original spin variable, the fluxes on the triangu- 
lar plaquettes correspond to the scalar spin chiralities 
S'l • S'2 X S'3, while (j)i2 is related to the 12 spins around 
the dodecagon loop. In the uniform spin liquid phase, 
5'L[0,0,0], fj^ = t{Oij,^ij) e R. For the non-uniform 
spin liquid state, the spin chirality arises, the fermion 
hopping will acquire a phase t^j-^ = t{6ij^(j)ij)e^^^^ ^ ren- 
dering a nonzero flux for a fermion moving around a loop 



E, 



loop • 



t(Oij^(l)ij)^ in principle, depends on the 



angles between spins Si and Sj. However in the mean 
field approximation, t{Oij^(j)ij) should be independent on 
the angles, ( 6>i j, ^^j) for the spin liquid phases because 
the fluxes through plaquettes are periodic in the whole 
lattice.^ Thus the hopping parameters t^j-^ can be classi- 
fied into two independent variables on the star lattice, tc 
labels the hopping amplitude inside the triangles and 
is the hopping amplitude between different triangles. For 



convenience, we introduce r = \ 
these two hopping amplitudes. 



to measure the ratio of 



where the effective hopping amplitudd^ 



cos(^)cos(^) + e" 



sm( — ) sm(- 



= t{Oij,(l)ij)e 



where the phase a^j is the vector potential generated by 
spin and corresponds to the Berry phase felt by the hop- 
ping electron. It is noted that the unit cell of the star 
lattice contains two triangular plaquettes and one 12-site 



For a given flux configuration, (j) = ^iQ^pdij^ there 
iaijd,m different phase configurations a^j. The effect of the 
phase configuration can shift the whole energy band in 
(3jthe k space, but do not modify the energy band structure. 
Namely, different choices of the phase configurations do 
not change the edge state and their topological properties 
of the star lattice. This allows us to set a simple choice 
of the phase configuration for various spin liquid phases 
to study their edge states and topological orders. 
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III. EDGE STATES AND TOPOLOGICAL 
ORDERS 



A. Bulk-edge conrrespondance theory of star 
lattice 



In general, we consider a strip of star lattice with the 
boundary along the ei direction and the periodic infinite 
62 direction shown in Fig.l. We assume that the spin 
liquid phase ^12], where the fiuxes satisfy the 

constraint, (/)<i + (/)i2 = 0. Using the Bloch theorem 



in the €2 direction, Cj 



j_ 

L2 



Y.j^e' '^Cniik), where n 



labels the the unit cell and I labels the sites in the unit 
cell. We set /c • 62 = k for convenience. The Hamiltonian 
in Eq.([2| can be written as 



Fe// = ^Ct(fc)Mfc)C(fc), 



(4) 



Where Ct(fc) = (c^(fc)...cl^6(fc)c^(fc)...4_6(fc)...4 ^(fc)), 
and 
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where Ni is the number of unit cell along the ei di- 
rection. The Bloch wave function can be written as 
— "iZn i^n,£cl^ £{k)\0) , where n runs the unit cell 
in the ei direction. Inserting it into the Schrodinger 
equation, = ^|^), the solution can be reduced to 

a set of equations (Harper equation) 



1,6 = ^V^n,l 



V^n,6e'^2 = ^^^^^ 
V^n,4 + + V^n,5e"*'^2 ^ ^^^^^ 



(8) 



where e = f-. Rewriting Eq.(^8| to a matrix form, we can 
express it in terms of a transier matrix form. 



(9) 



where M is a 2 x 2 matrix and its elements are 



Mn(£) 

Ml2(£) 

M2i(e) 
M22(e) 



• + -02 /a 
^ 2 

■ + -02 ^2 
-p* 2 

h 



where 



r^) cos ■ 



^2 

/a 
k 
k 



k 01 + ^2 , - 01 + 02 n 

2e cos - cos h cos 

2 2 2 ^ 

1 + r"^ - 2(1 + r^)e^ ^ - 2r^ cos k 

e{3 + 2r^ + r"^ - 2(2 + r^)^^ + e"^) - 2r'^e cos A: 



^3 + 2(1 

^3 + 2(1 



2r^ cos(/c + 0i) 
2r^ cos(/c — 02) 



We assume that the width Li of the star lattice contains 
an integer number of the unit cells, we can get the re- 
duced transfer matrix 



^Li + 1,1 



(M(e)) 



V^i,i 



(10) 



Considering the boundary condition '0Li,6 = '00 6 = 0, 
the edge energy Eedge satisfies {M{e)^^)2i = 0.1^ For 
Li ^ 1, the criterion for edge states follows^^ 

< 1 edge states localized in site 1 
I (M(/ij)) 11 1 ^ > 1 edge states localized in site Lx — 1 
1 coincide with bulk states 

(11) 

The quantum Hall conductance of systems can be ex- 
pressed in terms of the Chern number of U(l) bundle 
over the magnetic Brillouin zone.l^ The bulk-edge corre- 
spondence theory reveals that the Chern number C(/ij) is 
equivalent to the winding number of the edge state mov- 
ing around the hole of Riemann surface (the complex- 
energy surface), namely the intersection number between 
the canonical loop on the Riemann surface and the trace 
of the edge state energy /ij.^ ^ (see Fig. 3). Thus, the 
quantum Hall conductance can be given by the winding 
number of the edge states, cr^y^^ = —j^C{fij) when the 
Fermi energy lies in the jth energy gap. Therefore, we 
can count the Chern number from the energy spectrum 
of the system. 

The mean field studies of the star lattice with the 
Hamiltonian of Eq.Q give various spin liquid phases .^^ 
It is worth studying that the topological properties of 
these spin liquid phases. 
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FIG. 4: (Color online) The energy spectra of 6'L[0,0,0] for 
r = 1/2 in (a), r = 1 in (b), and r = 2 in (c). 



B. Uniform spin liquid phase: *S'L[0,0,0] 

The uniform spin liquid phase, [0,0,0] respects all 
the space-group symmetry of the lattice and time rever- 
sal symmetry. The energy spectra of the 6'L[0, 0, 0] phase 
for several set of parameters (r, (j)) are plotted in Fig. 4. 
It can be seen that both of the bulk energy band and edge 
states are k-symmetric, E{k) = E{—k)^ due to the space 
inverse symmetry, but the edge states are either embed- 
ded in the bulk states or isolated in the gap, namely there 
is no nontrivial bulk gap. Interestingly there exist two 
fiat bands lying in the energy-band gap and touching an 
edge- state band, which is caused by interference.^^ It is 
similar to a uniform spin liquid on the kagome lattice 
and can be spoiled by perturbations, such as the next 
nearest neighbor hopping.^ Thus, the Chern number is 
not well-defined and corresponds to common metals or 
insulators without IQHE. Actually the interactions from 
spinous can lead to instability of the uniform spin liq- 
uid phase to develop to the phases with breaking time 
reversal symmetry.^ 



C. Nematic spin liquid phase: 5'L[(/), — 0, 0] 

The spin liquid phases are characterized by a set of 
spin chirality operators. Different chiral spin phases ex- 
hibit different local magnetic fluxes. For the nematic spin 
liquid phase SL[<j)^ —<p,0]^ time reversal symmetry is bro- 
ken spontaneously,^^ we plot the energy spectra for some 
parameter (r, 0) in Fig. 5. It can be seen that the k- 
symmetries of both bulk bands and edge states are bro- 
ken, E{k) 7^ E{—k). It shows that time reversal symme- 
try breaking could induce space inverse symmetry break- 
ing that breaks the k-symmetries of bulk bands and edge 
states. However, the edge states are also either embedded 
in the bulk states or isolated in the gap. This implies the 
systems are common metals or insulators, but without 
IQHE, as [0,0,0] phase. 



D. Chiral spin liquid phase I: *S'L[0, 0, — 20] 

For the chiral spin liquid phase I 5'L[0, 0, — 20], time 
reversal symmetry is also broken spontaneously.^^ In 




(d) r = l^ = ^ (e) r = l,^ = - (f) r = h^ = ^ 

3 2 3 



FIG. 5: (Color online) The energy spectra of SL[—(j), 0, 0] for 
different parameters (r, 0) . 

principle, the Chern number in the jth. gap between the 
bulk energy bands depends on the parameters r and 0, 
Cj{r^(l)) (here we use this symbol for Chern number). 
However, we find from numerical investigations that the 
Chern number in the range of r and (j) obeys the following 
symmetries : 

(1) C^(r,27r- 0) = -C^(r,0) for G (0,7r); 

(2) Q(r,^ + (/)) = -C5-,(r,(/)) for G (0,^); 

(3) Q(r,(/)) = Q(-r,(/)) for G (0,2^); 

Thus, we can restrict the parameters only in the range 
of G (0, 1^) and r > 0. In Fig. 6 we plot the energy spec- 
trum for some typical parameters (t, (j)) and Li = 10. The 
Chern number can be counted by the winding number of 
the torus formed by two Riemann surfaces.^ (see Fig. 
3). The numbers in the right-hand side of each figures 
in Figs. 6 are the Chern number of the system when the 
Fermi energy lies in the corresponding energy gap. It can 
be seen that the bulk energy spectrum is k-symmetric, 
E{k) = E{—k) even though 7^ 0, but the edge state 
has no k-symmetry. It implies that the space inversion 
symmetry still holds due to the bulk band k-symmetry, 
but the spontaneous time reversal breaking breaks the 
edge state k-symmetry. Different Chern number implies 
different topological orders of the system. In order to 
give the phase diagram in the parameter space, we try to 
find out the critical lines in the parameter space. When 
the energy gaps close, the Chern number must change, 
namely a phase transition happens. Notice that the bulk 
spectrum has the Krammers degeneracy, E{k) = E{—k), 
it implies that the close of the energy gap of the bulk 
band happens only at the F point {k = 0). Thus, we 
solve the eigenenergies of the unit cell Hamiltonian (bulk 
band) at /c = 0. Numerical investigation indicates the 
bulk energy bands at /c = are linear with r, which al- 
lows us suppose that the eigenenergies of the bulk energy 
band at /c = have the form Ei{k = 0) = ±r + ^i(0), 
where i = 1,2,3. Substituting this form into the eigen 
equation of the Hamiltonian of the two-dimensional star 
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FIG. 6: (Color online) The energy spectra of 5'L[(/), 0, — 20] 
for different parameters (r, cj)) . 



lattice in the k space we can find that hi{(j)) satisfies a 
cubic equation, 



36 — 2 cos ( 







(12) 



When we consider G (0, f ), the solutions of the above 
cubic equation are 



hi =2cos^ 
63 = 2 cos 



2cos^ 



(13) 



which we set 61 < 62 < ^3- The closes of the bulk energy 
band gaps yield the critical lines in the parameter space, 



. 1 = V3si 



sm - 



sm 



\/3 
= >/3sin 



^ 3 



(14) 



These three functions separate the parameter space r — cj) 
into several different regions which have different topo- 
logical orders shown in Fig. 7. Different regions in the 
phase diagram represent the ground states with different 
Chern number configurations in which the Chern num- 
ber in different positions corresponds to the Fermi en- 
ergy in different gaps like Fig. 5. The phases with differ- 
ent Chern numbers reflect different integer quantum Hall 
conductance of the system. 



E. Chiral spin liquid phase II: 5'L[(/)i, 02, — (0i + ^2)] 

For more general chiral spin liquid phases II 
SL[(j)i^ 02, — ((/>i+(/>2)], the time reversal and space inverse 
symmetries are broken spontaneously. For example, the 
energy spectra of the spin liquid phase — 0] for 

some specific parameters are shown in Fig. 8. It can be 
seen that the k-symmetries of both bulk bands and edge 
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1.5 2.0 2.5 3.0 



FIG. 7: The phase diagram of 5'L[0, 0, —20]. Different region 
represents different Chern numbers for various filling fraction, 
which are even functions of r and periodic functions of with 
periodicity 27r. The energy spectrum along line 1 and line 2 
are shown along (a)^(d)^(b)^(c) and (d)^(e)^(f) in Fig. 
6 respectively. 



state are broken E{k) 7^ E{—k)^ and the Chern numbers 
in some energy gaps are not well-defined, such as the top 
three gaps for r = 1, 0i = ^ and 02 = Because the 
degenerate points of the bulk band and edge state are 
not at the k = point, the phase transition lines can not 
be solved easily. 



TABLE I: The energy band properties of different spin liquid 
phases 





EKS BKS TRS SIS 


Chern number 




r: 


1 
2 


1 2 


SL[0,0,0] 


yes yes yes yes 


-1 


-1 X 


■SL[f,-f,0] 


no no no no 


X 


X 


*^^l3 ' 3 ' ~"3~J 


no yes no yes 


-1 


-1 




no no no no 





-1 



EKS: Edge-state k-symmetry; BKS: Bulk band k-symmetry; 
TRS: Time reversal symmetry; SIS: Space inverse symmetry; 
x: non- well- defined Chern number. 
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FIG. 8: (Color online) The energy spectra of 
SL[(j)i, 02, —(01 + 02)] for different parameters (r, 0). 



IV. DISCUSSIONS 

To compare the basic energy-band properties of differ- 
ent spin liquid phases, we assume that the Fermi energy 
lies in the middle of the middle energy-band gap. The 
energy band symmetry and Chern numbers for some spe- 
cific cases are listed in Table I. 

It can be seen from Table I that for SL[0^ 0, 0] both of 
the edge k-symmetry (EKS) and bulk k-symmetry (BKS) 
are held due to the time reversal invariance and space 
inversion invariance. When r > 2, the Chern number 
becomes non-well-defined. The cases in the 2 ~ 4 lines in 
Table I indicate that there is no EKS, but BKS remains 
for 1^, — This implies that the spontaneous 

time reversal symmetry breaking does not always induce 
the space inversion symmetry breaking and the EKS can 
be broken by the time reversal symmetry. 

The results in Table I reveal that the Chern num- 
ber depends on not only the time reversal and space in- 
verse symmetries, but also the parameters (r, 0i, 02, ^12) 
of the star lattice. The ground states become nor- 
mal metals or semiconductors for the phases without 
the well-defined Chern number. Interestingly, there is 
a topological invariance for the exchange of the mag- 
netic fluxes in the two triangles 5'L[0i, 02, — (0i + 02)] 
and S'L[02,0i, —(^1 + 02)]- They have the same Chern 
numbers and their bulk energy bands are k-asymmetric. 
These findings indicate some new phases in the 2- 



dimensional materials.!^ 



V. CONCLUSIONS 

In summary, we have studied the edge states and their 
topological orders in the different spin liquid phases of 
star lattice by using the bulk-edge correspondence theory. 
The bulk and edge-state energy structures and Chern 
number depend on the spin liquid phases and hopping 
parameters because the local spontaneous magnetic fiux 
in the spin liquid phases breaks the time reversal and 
space inversion symmetries. We have given the charac- 
teristics of bulk and edge energy structures and their cor- 
responding Chern numbers in the uniform, nematic and 
chiral spin liquid phases. In particular, we have obtained 
analytically the phase transition lines of different topo- 
logical phases and their corresponding phase diagrams for 
the chiral spin liquid states SL[(j)^ 0, —20]. We have also 
found that the topological invariance for the spin liquid 
phases,6'L[0i, 02, -(^1 + ^2)] and 5'L[02, 0i, -(^1 +^2)]. 
The results tell us the relationship between the energy- 
band and edge-state structures and their topological or- 
ders of the star lattice. Especially, this star lattice 
has been synthesized in the material called iron acetate 
recently. Therefore, our results provide a Hall conduc- 
tance experimental guideline to discriminate the spin liq- 
uid phases in real materials and cold atoms in optical lat- 
tice. The changes of filling fraction could be implemented 
by tuning the applied gate voltage. These results can also 
give some hints for understanding the Heisenberg model 
on the star lattice. 
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